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The aim of this paper is to derive an analytical expression for the self-propulsion veloc- 
ity of a micro-swimmer that consists of N spheres moving along a fixed line. The spheres 
are linked to each other by the rods of the prescribed lengths changing periodically. 
For the derivation we use the asymptotic procedure containing the two-timing method 
and a distinguished limit. Our final formula shows that in the main approximation the 
self-propulsion velocity is determined by the interactions between all possible triplets of 
spheres. 



1. Introduction and formulation of problem 

1.1. Introduction 

The studies of simple micro-swimmers (or micro-robots) represent a flourishing modern 



research topic (see Purcell (1977) Alexander, Pooley and Yeomans (2009) Earl, Pooley, Ryder, Bredberg and Yeomans 



Alouges, DeSimone and Lefebvre (2008)| Becker, Koelher, and Stone (2003) Gilbert, Ogrin, Petrov, and Winlove (2011 



Golestanian & Ajdari (2008) Golestanian & Ajdari (2009) I, which creates the funda- 



mental base for modern applications in medicine and other areas. In this paper we gen- 



eralize the three-sphere micro-swimmer by Golestanian & Ajdari (2008) to the general 
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TV-sphere micro-swimmer. It is possible since the employed two-timing method and a 
distinguished limit significantly simplifies the analytical procedure. 

1.2. Formulation of problem 

We consider a micro-swimmer consisting of N rigid spheres of radii R*, i = 1, 2, . . . N 
with their centers at the points x*(t*) of the x*-axis (x* +1 > x*); t* is time. The spheres 
are connected by TV — 1 rods of lengths if* = x*. — x* where our choice is always k > i. 
The masses of the spheres and the rods (in the Stokes approximation) are zero. The rods 
are so thin that their interaction with a fluid is negligible. The lengths of the rods are 
prescribed as 

Zf=£f+Af (1.1) 

where L\* are the averaged values and Xf*(r) are the oscillations, which are prescribed as 
27r-pcriodic functions of r = tot* with a constant frequency lo. Asterisks mark dimensional 
variables and parameters. 

In the Stokes approximation the total force acting on each sphere is zero (their masses 
are zero) , hence the equation of motion for the i-th sphere can be written as 

<z*-^3<i?^/(2Zf) = -/* (1.2) 

where k* = 6irr)R*, rj is viscosity, dots above the functions stands for d/dt*. The 
l.h.s. of (|1.2[) represents a viscous friction, while /* is the force exerted by the rods 



to the i-th sphere. In order to derive (|1.2j) we have used the fact (see Lamb (1932) 



Landau & Lifshitz (1959) Moffatt (1996) I that a sphere of radius R*. and position x*. 
moving along the x*-axis with velocity x* k creates at the center of «-th sphere the x*- 
component of fluid velocity equal to —3R*,x*./ (21!?*), where the minus sign corresponds to 
x± > Xj. The considered mechanical system is a closed one, hence the total force exerted 
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N 



E/* = ° (1-3) 



Eqns. (|1.2[) , (|1.3p represent the system of ODEs to be solved in this paper. Notice that 
we do not use the summation convention. 

The equation (|1.2[) and its solution x*(t*) = {x\, x\, . . . , x^), contain three character- 
istic lengths: the radius R of spheres, the distance L between the neighbouring spheres, 
and the amplitude A of oscillations of rod lengths; at the same time the only explicit 
characteristic time-scale T corresponds to the frequency ui: 

R, L, A, T=l/u (1.4) 
The dimensionless variables and small parameters are 

x* = Lx, L* k = LL\, R* = i?i?. t , Af = AAf , t* = Tt (1.5) 
/* = -QirriRLfr/T, e = X/L < 1, 5 = 3R/(2L) < 1 

Then the dimensionless ecms. (|l.l[) - ([1.3p take the form 

Ri±i — 5 ^ Rikik/hk = fi, Rik = RiRk (1-6) 

hk = Lik + eXik (1.7) 

J2fi = f-I = 0, I=(l,l,...,l) (1.8) 

i 

One should note that 'dots' above function in (|1.2[) and (| 1 .€>[) correspond to the dimen- 
sional and dimensionless time derivatives correspondingly. The first equation (|1.6[) can 
be rewritten in the matrix form 

N 

Ax = f or ^2 A ik£k = fi (1-9) 



k=l 




for i = k, 

(1.10) 



-SRik/hk for i ^ k 
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1.3. Notations 

The variables x — (x\,X2, ...,Xn), t, s, and r serve as dimensionlcss coordinates of 
spheres, physical time, slow time, and fast time. We use the following definitions and 
notations: 

(i) A dimensionless function / = f(s, r) belongs to the class 0(1) if / = 0(1) and 
all partial s-, and r-derivatives of / (required for our consideration) are also 0(1). In 
this paper all small parameters appear as explicit multipliers, while all functions always 
belong to D(l)-class. 

(ii) We consider only periodically oscillating in r functions 

feSi: f(s,r) = f(s,T + 2n) (1.11) 

where the s-dependence is not specified. Hence / G Sj f)D(l). 

(iii) The subscripts r and s denote the related partial derivatives. 

(iv) For an arbitrary / £ $j the averaging operation is 

1 /.TO+27T 

</> = 2^y /M^> Vt "o (1-12) 

where during the integration we keep s = const and (/} does not depend on r . 

(v) The class of tilde-functions (or purely oscillating functions) is such that 

/: 7( S) r) = /( S ,r + 27r), with (/ ) = 0. (1.13) 

Tilde-functions represent a special case of i^-functions with zero average. 

(vi) The class of bar-functions (or mean-functions) is defined as 

/: / r ^0, 7( S ) = (7( S )) (1.14) 
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2. Asymptotic procedure 

The use of e-dependence of (jl.ip leads to the presentation of matrix A (jl.lOj) as a 

series for e — >• (we consider 5 as a fixed parameter) 

A = C + eSA' Q + C =A + («o (2.1) 

for i = k, 

RikKk/L 2 lk fori^k 
where we do not present the expression for Bo since it will not affect the answer. 

The crucial step of our procedure is the introduction of a fast time variable r and a 
slow time variable s. We take r = t (which corresponds to the prescribed oscillations of 
the rods) and s = e 2 t. This choice of s can be justified by the same distinguished limit 
arguments as in Vladimirov (2012)] here we present this fact without proof, referring 



k = diag{i?i,i? 2 , ...,R N }, 



only to the most important fact that it leads to a valid asymptotic procedure. Therefore 
we use the chain rule d/dt = d/dr + e 2 d/ds and then we accept (temporarily) that r 
and s represent two independent variables. The two-timing form of eqn. (|1.9[) is 



(C +eSA' + ...)(x T +e 2 x s ) = f (2.2) 

where unknown functions are taken as the series 

x(t,s) =x (s)+sxi(t,s) + f(r,s) = f (r, s) + ef^r, s) + . . . (2.3) 

The accepted condition Xq = reflects the fact that the large distances of self-swimming 
are driven by small self-oscillations. Now we consider the successive approximations of 

(E2D,<E31) in e: 

(i) Terms 0(e°) give f = 0. 

(ii) Terms 0(e 1 ) give Coa;or = fi', the averaged part of this equation gives f 1 = 0, while 
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the oscillating part yields 
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'0 x 1t — fl 



(2.4) 



(iii) Terms 0(e 2 ) give the equation CqX2 T + SA' Q Xi T + CqXq s = / 2 ; its averaged part is 



l x 0s + 5(A' x lT ) = f 2 



The force f 2 can be excluded from (|2.5p by (|1.8I) : 



I • Co^ 0s + SI ■ (Kx lT ) = 



(2.5) 



(2.6) 



The self-propulsion with averaged velocity Vq means that Xq s — VqI, hence 

Vo = -f J -j^> (2.7) 



I-AqI 

where the matrix Co is replaced with Ao since we consider only the main (linear) term in 
5. Expression (|2.7p still contains unknown functions X\ T which can be determined from 
(|2.4[) with the use of constrains (|1.7[) . (|1.8[) . Indeed, the equation (|2.4[) (with linear in S 
precision) gives X\ T = Ap f-y with A^ 1 = diagjl/i?!, l/i?2: ■ • ■ , ~^/Rn}\ it means that 
Xir — g with the components cji = fu/Ri. One can see that (jl.7|) yields ^ — Iji = \f T , 
while (|1.8|) leads to ^ i?i<7i = 0. Both these restrictions can be written as one N x N 
matrix equation 



Mg = l 7 



-1 1 

-1 1 

0-0 

i?l i?2 R3 



-1 1 

Rn-i Rn 



I = 



4 



A AT-1 



V J 



(2.8) 



The substitution of x lT = M 1 l T into (|2.7p gives us self-propulsion velocity in the matrix 
form 



Vo = -5- 



(2.9) 
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where the matrix 



is 



(_l)^+i AM -i = 



Ai — A 


A 2 - 


A 


A, 


-A .. 


■ AjV-1 


-A 


1 


Ai 


A 2 - 


A 


A 3 


-A .. 


• AjV-1 


-A 


1 


A a 


A 2 




A, 


-A .. 


■ AjV-1 


-A 


1 



Ai 
Ai 



A 2 
A 2 



A 3 
A 3 



Ajv-i-A 1 



A 



N-l 



(2.10) 



1 ) 



A fc = J2 R ^ k>l; A = A N . 



Further calculations show that (|2 can be presented as 

V ° = A2 Gikl 



(2.11) 



i<k<l 

dki = RiRkRi 



T 2 T 2 I {^ik^klr — Aife T Afej) 



1 1 1 

T ' 7~2 72" 

ik n kl u il , 

where the sum is taken over all possible triplets (i, k, I) : l<i<k<l<N. 
For the three-swimmer this sum contains the only term, which coincides with one by 



Golestanian & Ajdari (2008) In general it contains N\/[(N — 3)!3!] terms: for the four- 
swimmer we already have four triplets (1, 2, 3), (1, 2, 4), (1, 3, 4), (2, 3, 4), for the five- 
swimmer - 10 terms, while for the ten-swimmer the number of triplets grows up to 
120. 

The expressions for M^ 1 (|2 . 10[) and Vq (|2.11[) have been obtained by the explicit 
calculations for N = 3, 4, 5 and by the mathematical induction for any N. 

Formula (|2.11[) represents the main result of this paper. According to (|2.1ip I ■ x s = 
Vq = 0(5); however physical velocity is I ■ Xt = s 2 Vq. Hence the order of magnitude of 
the dimcnsionlcss physical velocity is 0(e 2 6). 
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3. Discussion 

1. The explicit formula (|2.11j) allows one to find the optimal strokes, to calculate the 
required power, the efficiency of self-swimming, and all related forces (both oscillatory 
and averaged). However the large number of terms in (|2.1ip makes all these problems 
rather cumbersome, and places them out of the scope of this short paper. 

2. Our approach (based on the two-timing method and distinguished limit) is techni- 
cally different from all previous studies in this area. The results for iV-sphere swimmer 
show its analytical strength. 

3. The expression (|2.11l) can be predicted without any calculations, on the base of 
the result for TV = 3. Indeed, if we are interested in the main term of the order e 2 5, 
then only the triple interactions can be taken into account, as they have been described 



by Golestanian & Ajdari (2008) The additional (to triplets) interactions between four 
spheres will inevitably produce the next order term 0(e 3 8), which we do not consider. 

4. There are some interesting discussions about the physical mechanism of self-propulsion 
in the quoted literature. However one can also notice that a similar result does exist for 



self-propulsion in an inviscid fluid (Saffman (1967) I and some physical explanation can 



be achieved if we replace the term 'virtual mass of a dumbbell' by the term 'viscous 
drag coefficient of a dumbbell'. Say, for a three-sphere swimmer this coefficient decreases 
when the distance between two neighbouring spheres (a dumbbell) decreases and then 
the third sphere is used to 'push' or 'pull'. If the reverse motion of the third sphere meets 
the increased drag coefficient of the dumbbell, then self-propulsion is achieved. 

5. The mathematical justification of the presented results by the estimation of the error 



in the original equation can be performed similar to Vladimirov (2010) Vladimirov (2011) 



6. One can also derive the higher approximations of self-propulsion velocity, as it has 



been done by Vladimirov (2010) Vladimirov (2011) They can be especially useful for 
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the studies of motion with Vq = (say, if all correlations involved to (|2.11[) arc zero). In 

this case one can show that self-propulsion can be generated by interactions of four and 

more spheres. 

7. In this paper we consider only periodic oscillations of constraints. The studies of non- 
periodic oscillations might represent an interesting problem. An attempt in this direction 



have been made by Golestanian & Ajdari (2009) In fact, such generalizations have been 



already considered for many different oscillating systems. 

The author is grateful to Profs. R. Golestanian, A. D. Gilbert, and H.K.Moffatt for useful 
discussions. 
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